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It was shown recently that the anomalous scaling of simultaneous correlation functions in tur- 
bulence is intimately related to the breaking of temporal scale invariance, which is equivalent to 
the appearance of infinitely many times scales in the time dependence of time-correlation functions. 
In this paper we derive a continued fraction representation of turbulent time correlation functions 
which is exact and in which the multiplicity of time scales is explicit. We demonstrate that this form 
yields precisely the same scaling laws for time derivatives and time integrals as the "multi-fractal" 
representation that was used before. Truncating the continued fraction representation yields the 
"best" estimates of time correlation functions if the given information is limited to the scaling ex- 
ponents of the simultaneous correlation functions up to a certain, finite order. It is worth noting 
that the derivation of a continued fraction representation obtained here for an operator which is not 
Hermitian or anti-Hermitian may be of independent interest. 



I. INTRODUCTION 

It is commonly argued that fully developed hydro- 
dynamic turbulence exhibits simultaneous statistical ob- 
jects whose scaling properties are anomalous. For exam- 
ple the so called structure functions satisfy scaling laws 
of the form 

SniR) = {\uir + R,t)-u{r,t)n^R'^- , (1) 

where u{r,t) is the Eulerian velocity field, and are 
scaling exponents which are nonlinear functions of n. The 
nonlinear dependence is referred to as "anomalous scal- 
ing" or "multi-scaling" , and the issue of evaluating these 
exponents from either phenomenological models or from 
first principles has attracted significant amount of effort 
in the last decade. 

It has only recently been discovered that also the 
time dependence of the n-th order correlation functions 
is multi-scaling, and that "dynamical scaling" is broken. 
This phenomenon seems to distinguish turbulence from 
other problems in which scaling is anomalous, like crit- 
ical phenomena. In the latter case dynamical scaling is 
invoked by stating that a space time correlation function 
F{R, t) is a homogeneous function of its arguments in the 
sense that F{XR, XH) = \'^F{R, t), where C and z are the 
"static" and "dynamic" scaling exponents respectively. 
In turbulence such relations do not exist even when the 
same-time correlation functions are homogeneous func- 
tions of the spatial coordinates. The importance of this 
fact in determining the structure of the theory has been 
stressed in [Q, and see also [1). 

In this paper we address temporal multi-scaling on the 
basis of the continued fraction representation of turbulent 
correlation function |6|j^ . This approach will lead us to a 
different point of view of temporal multi-scaling, in agree- 
ment with the conclusions as Ref. Q. The advantage of 
the present formulation is three-fold: first, it is derived 
on the basis of an exact formulation of the time correla- 
tion functions and their dynamics. The phenomenon of 



temporal multi-scaling is related in this approach to the 
scaling properties of higher order temporal derivatives of 
correlation functions, computed at zero time. Second, 
this approach furnishes information not only about lead- 
ing scaling exponents, but also about sub-leading ones. 
Third, a finite truncation of the continued fraction rep- 
resentation is in a sense the "best" possible representa- 
tion when the information about the scaling of one-time 
correlation functions is limited to the low order scaling 
exponents. We will show that the scaling laws exhibited 
by the continued fraction representation are identical to 
those predicted by the "multi-fractal" representation of 
time-correlation function |^, adding weight and justifica- 
tion to the latter. Since the multi-fractal representation 
was used recently to estimate the scaling exponents from 
first principles , we ascribe some weight to being able 
to justify it further. 

To keep the formalism minimal and the result clearest, 
we treat in this paper only the second order space-time 
correlation function of turbulent fields. The formalism 
can be used to generate representations of any higher 
order correlation function, but we do not elaborate on 
this in the present text. In Section II we review briefly 
the Zwanzig-Mori formalism [^|j9| which was applied to 
time correlation functions in turbulence ^,|^ , and display 
the continued fraction representation of the second order 
time correlation function of the velocity field. We show 
that the coefficients in this representation can be writ- 
ten in terms of time derivatives at zero time of the same 
second order correlation function. In Section III we set 
up a theory for the evaluation of the scaling exponents 
of these coefficients. In Section IV we derive the scaling 
laws implied by the continued fraction representation for 
the time derivatives of the correlation function evaluated 
at time zero. In Section V we relate our results to the 
multi-fractal representation of the correlations functions, 
and explain the scaling-equivalence of the two represen- 
tations. In Section VI we offer a summary and a short 
discussion. 
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II. CONTINUED FRACTION REPRESENTATION 



where 



In thinking about dynamics one cannot deal with time- 
correlation functions of the Eulerian field since these are 
dominated by the kinematic sweeping time scale. We 
need to consider Lagrangian or Belinicher-L'vov velocity 
differences. We prefer the latter since they obey Navier- 
Stokes like equations of motion which are local in time. 
In terms of the Eulerian velocity u{r, t) Belinicher and 
L'vov defined ||l^] the field v{rQ,to\r,t) as 



v{ro,to\r,t) 



u r - 



■p{ro,t),t] 



(2) 



p(ro,i)= / dsu[ro + p(ro,s),s] 

J to 



(3) 



The observation of Belinicher and L'vov was that the 
variables W{rQ,tQ\r,r,t) defined as 

W(ro, to|r, r',t)= v{ro,to\r, t) - t;(ro, io|r', t) , (4) 

exactly satisfy a Navier-Stokes like equation in the in- 
compressible limit: 



dt 



- PW(ro, h\r, ro, t) ■ V, + P'W(ro, to|r', ro, t) ■ - + V;^) W(ro, iok, r', t) = 



(5) 



We remind the reader that the application of the 
transversal projector V to any given vector field a(r) 
is non local, and has the form: 

[Pa(r)]" j drV^^ir - f)a^(f). (6) 

The explicit form of the kernel P°'^{r) can be found, for 
example, in [Q. In (||) P and V are projection opera- 
tors which act on fields that depend on the corresponding 
coordinates r and r' . For our purpose below we intro- 
duce the Liouville operator C{ro, to\r, r') which is defined 
via the total time derivative of W(ro, to\r,r' ,t) at time 



dW{ro,to\r,r',t) 



dt 



C{ro,to\r,r')yViro,to\r,r' ,to) ■ 



(7) 



The identity of the Liouville operator follows by defini- 
tion. 



dv{ro,to\r - p{rQ,t),t\,t) _ du{r,t) 



du{r, t) 



dt 

+ [u{r,t) ■ V]u{r,t) 



dt 



(8) 



Translating all coordinates by p(ro,t) we find after ele- 
mentary algebra. 



d 

C{ro,to\r,r') = — 



+ v{ro,to\r,tQ) ■ V,. 
+ v{ro,to\r' ,to) ■ V,.- 



(9) 



Consider now the time dependence of "fully fused" 2nd 
order correlation function 

Tfiro\r, r', r) = (W"(ro|r, r', to)W''iro\r, r', to + r)) . 

(10) 



By "fully fused" we mean here that the coordinates of 
the two velocity differences are the same, and they differ 
only in their time argument. The same-time counterpart 
of this correlation function, i.e. ^2{'''o\r,r',T = 0) is 
independent of ro and in an isotropic and homogeneous 
ensemble it is a function of |r — r'| only. Accordingly it 
differs from the standard structure function 5*2 (|t — r'\) 
only in having the full 2nd rank tensorial character. For 
this analysis we choose all the three vector distances to 
have the modulus in the inertial range, of the order of R. 
In the notation we will keep only this R but remember 
that the angular dependence is not shown explicitly. 

To proceed, we consider the tensorial correlation func- 
tion ( px| ) as an inner product in the space of vectors W, 
denoted as 



(11) 



In particular we are interested in the Laplace transform 
ofEq. (0) 



T2{R, z) = iW, 



1 



:W) 



(12) 



It has been shown by Grossmann and Thomae M that 
the Zwanzig-Mori projection operator formalism MM ap- 
plies to turbulent systems described by Navier-Stokcs like 
equations. In it has been demonstrated that the con- 
tribution of the memory kernel and higher order con- 
tinued fraction is considerable. The central idea is to 
decompose the resolvent 



7^(2) 



1 



c 



(13) 



by means of a projection operator P. With Q = 1 — P as 
projector orthogonal to P one has the resolvent identity 



2 



pn{z)p = p- 



1 



PCP - PCQ 



1 



z-QCQ 



QCP 



(14) 



Since we note that the time correlation function of inter- 
est is the W — W matrix element of the resolvent, it is 
useful to choose P to be the projector on W, 



P. = (w,.)(w,yv)"^w 



(15) 



Its basic properties are PP = P (idempotent) and 
P^ = P (self adjoint), characterizing an orthogonal pro- 
jection. ^From Eq.(|l^ we have an expression for P: 



p = pn{z)p 



PCP - PCQ- 



QCQ 



QCP 



(16) 



Compute now the W — W matrix element, use PW = 
W, as well as the definition of P to compute 

(w,7^(z)w)(w,w)^M^(w,w) (17) 

-(W,£W) - {W,CQ i-j^QCW)] - (W,>V) . 

z - QCQ 

We therefore conclude with 

fco(i?) 



:F2iR,z) = 



z-j„{R)-Ko{R,z) 



(18) 



where 



fco(i?) = (W,W), 

7o(i?) = (W,^W)(W,W)-S 

ko{R, z) = (Qt£tw, L^QCW) (W, W) 

^ z - QCQ ' 



(19) 



Here, of course, = Q, but in the next steps of gen- 
erating the continued fraction hermiticity will not hold. 
Realizing that the kernel K.^ {R, z) has the same resolvent 
structure as J^2{R,z) except that it now features QCQ 
instead of C and the states are Q^OW, QCW instead of 
W, W, with C^ the adjoint of £ - one can continue the 
fraction by the same procedure. This is more transparent 
if we denote 



Wi = Qcw, yvi = Q^c^yv, Ci = qcq 

so that Ko{R, z) takes on the form 



ko{R,z) = Wi, 



1 



Wi (W,W)" 



z — Ci 

Now we define a new projection operator 

A - = (Wi,-)(Wi,Wi)-i>Vi . 



(20) 



(21) 



(22) 



As a result of Wi being different from Wi when C is 
not Hermitian or anti-Hermitian, this operator Pi is not 



Hermitian and performs accordingly non-orthogonal pro- 
jections. But Pi still is idempotent, PiA = A) which 
is the essential property for deriving the analogous resol- 
vent identity with Ci in (|2l| ) as for the original resolvent 
(|l^ ) with C. Defining Qi = 1 — A we can repeat the 
argument leading to (O) and dl^), and find 



Ko(i?,z) 



ki{R) 



z~-ii{R)-Ki{R,z) 
where 

ki{R) = (Wi,Wi)(W,W)-S 
7i(i?,) = (VVi,AWi)(VVi, Wi)-S 



(23) 



(24) 



1 



Ki{R,z)^{Q\C\yVi, 

^ z - QiCiQi 

x(>Vi,Wi)-i . 

Hence starting from 

1 



QiCiW, 



k,,{R, z) (>v„+i, — — yy„+i ) (vv„, w„)-i 

^ Z - Cn+i ' 



with n > one arrives at 



kn{R,z) 



kn+l{R) 



(25) 



(26) 



Z - 7n+l (R) - Kn+l{R, Z) 

where 

fe„+i(i?) = (VV„+i,>V„+i)(W„,W„)-i , 

7„+i(^) = (VV„+i,£„+iW„+i)(W„+i,W„+i)-i . (27) 

Here we used the notation 

W„+l = QnCnWn , VV„+1 = QI^Wu , 

£„+i = QnCnQn , Wo - Wo - W , £o - ^ , (28) 

and defined new projection operators as 

A • = (vv„,-)(vv„,w„)-iw„ 

Q« = i - A . (29) 

From (|l^) and ( p6| ) it thus follows that the Laplace trans- 
form J^2(R, z) of the correlation function ( |l0| ) can be 
written in continued fraction representation: 

ko{R) 



:f2{r,z) = 



z--foiR) 



kiiR) 



7i(A- 



fc2(A 



72 (^) 



(30) 



The novelty when the operator C is not Hermitian or 
anti-Hermitian as is the case here is that the new projec- 
tion operators introduced to continue the fraction per- 
form non-orthogonal projections. Only the lowest order 
iterates of A and Qn, i.e. the operators P and Q, per- 
form orthogonal projections. In the following section we 
analyze the scaling properties of this representation. 
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III. SCALING PROPERTIES OF THE 
CONTINUED FRACTION REPRESENTATION 

In this section we determine the leading scahng expo- 
nents of the coefficients which appear in the continued 
fraction representation (^0|). This leading scahng be- 
havior is given in terms of correlation functions of time 
derivatives of the velocity differences W computed at 
zero time: 

fc„(i?) « (£">v£"w)/(i:"-i>Vi:"-^w) , (3i) 
7„(i?) « (£"+iWi:"W)/(i:"W£"w) . 

Here the symbol « means "leading scaling order". 
Eq.(^l[) yields the following explicit scaling 



fco oc R'^^ , 

7„(i?) CX i?C2,. + 3-C2„ + 2-l ^ 



(32) 



The rest of this section is a demonstration of this re- 
sult. The reader who prefers to see the connection to the 
multi-fractal representation can go directly to the next 
section. 

Let us define 



J n 



■ Qn-lQnQn-1 ■ ■ ■ Qo 



Now 



(33) 
(34) 

(35) 



implies that 

Jn = QoQl ■ ■ ■ Qn = QnQn-l ■ ■ ■ Qo ~ Qn 
^l-Po-...-Pn ■ 

Hence Eq. (|8|) leads to 

Wn+l = J„CJn-lCJn-2 ■ ■ ■ CJqCVS; , 
Cn+iyVn+1 = JnCJnCJn-l ■ ■ ■ CJqCW . (36) 

It follows that 

(VV„+i,W„+i) 

= (vw, cjoCJi ■ ■ ■ cjn-iCJnCJn-i ■ ■ . cjotyv) , 

(W„+i,£„+iW„+i) 

= (W, LJqCJi ■ ■ ■ CJnCJnCJn-1 ■ ■ ■ CJqCW) . (37) 

In order to arrive at Eq. (^Tj) we have to show that 

(VV„+i,W„+i) « (W,£2"+2W) 

(VV„+i, £„+iW„+i) « (W, . (38) 

It suffices to use the following assertion which we prove 
below 



(>V,£'=Jg/gW) « (W,£*'7,W), k>q>Q, (39) 

where Iq is any composition of C and Ji which has two 
properties. First, it contains at least q + 1 times the op- 
erator C which we indicate by the subscript q in Iq and 
second, on the right hand of each Ji there are at least 
i + 1 operators C in the composition. Since the compo- 
sitions of C and Ji appearing in ( ^7|) satisfy these two 
properties, (^8|) results straightforwardly. 

We now proceed to the proof of ( |39| ) by induction. As- 
suming for p < q < k that 



iw^c' jpipW) ~ (w,£'=/pW) , 

we show that 



(40) 



(W, C'^JqiqW) « (W, C'^iqW) . (41) 

Using Eq. (||) the LHS of Eq. (||) becomes 

{W,l''JqIqW) = {W,l''Jq-liqW) - {W , I'' PqlqW) . 



(42) 



By Eq. ( [40D the first term on the RHS scales as 

(W, C^'Jq-liqW) - (W, C'^iqW) . (43) 

As for the second term, without specifying Iq we notice 
using Eq. ( |29| ) for Pq that it is the product of the fol- 
lowing three factors that feature only operators Jp with 
p < q and with more than p times the operator C on their 
right hand so that by (40) one has 



(W, C''Wq) = (W, tjq-lCJq-2 ■ ■ ■ CJoCW) 

« (W,£'=+«W), 

{Wq, Wq) = (W, tjoljl . . . CJq-l . . . CJoCW) 

« (W,>C2«W), 

iWqJqW) - (W, CJoCJl . . . CJq^liqW) 

« (W, tnqW) . (44) 

In order to compare the scaling behavior of the two 
terms on the RHS of Eq. (42) in the limit of infinite 
Reynolds number we recall that within the inner prod- 
uct (or equivalently the average operation) an operator C 
amounts to simply W/i? (3j. This follows from the con- 
vergence in the UV and in the IR (in the limit Re— > oo) 
of the integral implied by the terms VW ■ in Eq. (||), 
so that the leading contribution comes from distances of 
the order of R. In other words, £ in a correlation func- 
tion, when it operates on W, introduces a term of the 
order of W • VW, which, due to the demonstrated lo- 
cality in scale space, can be estimated as adding to the 
correlation a factor of the order of W/i?. Hence one has 
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(45) 



where the last step follows from Eq. (|^). 

The composition Iq containing say j* times the opera- 
tor C with j* > q and having at least i times the operator 
C on the right hand of each operator Ji, one can write 



3 



J=0 



(46) 



where Cj{R) is a function of the separation distance R 
resulting from all the contributions to OW. Note in 
particular that Cj-{R) = 1. Considering now the follow- 
ing quantity for arbitrary s one has 



(w,/:./,>v) cx 



{R)R^^ 



(47) 



Because of the Holder inequalities 

Cm - Cm-r < C" - Ci-r m> n , T > , (48) 

to which any two contributions of Eq. ( ^7[ ) can be re- 
duced, the leading contribution is that for which the in- 
dex of C is maximum, i.e. s + j* + 2. As a result 

(W, CsiqW) (X rCs+j-+2-i'-^ . (49) 

It follows from Eqs. (|4^) and ( ^ ) that the scaling be- 
havior of the first term on the RHS of Eq. (42) is 



(50) 



whereas using Eqs. (p4|),(p5[) and ( |4^ ) that of the second 
one is 



■k-j' 



(51) 



Recalling that k > q and j* > q one has by Eq. ( ^8| ) 

Cfc+j*+2 - Cfc+g+2 < Cq+j'+2 - C29+2 ■ (52) 

Hence expression (|o|) which coincides with the RHS of 
Eq. (41) is the leading one. We have therefore proven 
that relation ( pO| ) implies relation (|4l|). Now it remains 
to show that Eq. (^) is true for p = 0. But this follows 
straightforwardly by setting q = in Eq. (^) and re- 
placing J_i by 1. Hence assertion ( p9| ) is proven. This in 
turn implies as was shown at the beginning of this section 
that (Sll) is satisfied. 



IV. SCALING LAWS IMPLIED BY THE 
CONTINUED FRACTION REPRESENTATION: 
DERIVATIVES AT TIME ZERO 



To prepare for the comparison between the continued 
fraction and the multi-fractal representations we identify 
in this Section the leading scaling exponents that char- 
acterize the nth order time derivative of the correlation 
function Eq. ([ll|) at r = 0. We show that 



Here we use the shorthand notation 

d"J^2{R,T) 



^2 (R^o) 



(53) 



(54) 



r=0 



In the next Section we will show that the same exponents 
are predicted by the multi-fractal representation, making 
the prediction of the Taylor expansion of the correlation 
function the same from the point of view of scaling be- 
havior. 



From Eq. (18) one deduces by inverse Laplace trans- 



form the following equation for .^2^''(i?, r) 

T^^\r,t) = ^^{R)T2{R,t) (55) 

+ / Ka{R,T')T2{R,T-T')dT' , 

Jo 

where Kq{R,t) is the inverse Laplace transform of 
ko{R,z) [Eq. (IH)], 



Ko{R,t) 



(>Vi,e^i"Wi) 



fco(i?) 



(56) 



Equation ( [Sq ) is the so-called memory-function equation, 
Kq{R,t) being the memory kernel. At r = Eq. ( [55| ) 
becomes 



J^i^\R,Q) = Jo{R)ko{R) cx R^^'^ 



(57) 



where we used the scaling laws (P2|). 

The higher order partial time derivatives are obtained 
by differentiating Eq. (pa), 



T^^\R,r)^j,iR)Ti--'\R,r) + Y,Ki-'-'\R,r)Ti'\R,0)+ / dT'Ko{R,T')J-^r'\R,r-r') . (58) 



5 



At r = one obtains 

^("'(i?,0)=7o(i?)^^""'Hi?,0) 

n-2 
k=0 



(59) 



Now one realizes that this leads us to consider a hierarchy 
of equations for K'f\R,0), q > 1, i > 0- Noticing that 
Eq. ( p^ ) is formall y t he same as Eq. (|8|) one obtains in 
analogy with Eq. ( p5|) 

Kf^(i?,T) =7,+i (i?)K,(i?,T) (60) 

+ / K,+ l{R,T')K,{R,T-T')dT' , 
JQ 

where Ki {R, r) is the inverse Laplace transform of Eq. 

K,{R,t) = (W,+i,e^'+i^W.+i)(W„W.)"' . (61) 
Notice that at r = using also Eq. (|2^) yields 

K,{R,0)^k,+i{R) . (62) 
In analogy with Eq. ( |59| ) one has for q > 1 

K['>\R,0)^j,^,iR)K['^-'\R,0) (63) 

+ k[1-'-'\r,o)k\'\r,o) . 

k=0 

Note that deno ting from now on K i = ^2 and allowing 
j > — 1 in Eq. (63) one recovers Eq. (|59|). 

It follows that to determine the scaling behavior of 
k'i[^^ (_R, 0) for n > 1 in order to prove Eq. ( ^ ) one has 

to know K'i''^ {R,0) V q + 2i < n - 2 with g > and 
i > —1. It is sufficent to prove the following 



(64) 



where g > 0, i > —1 and Si^-i is a Kronecker delta, 
which we do by induction in the sequel. Setting i = — 1 
and recalling that Co = Eq. (|5^) results directly. 

Let us assume that Eq. (|6^) is true Vq + 2j<n — 3 
with n > 1. Then we show that this expression is still 
valid for n + 1. By the scaling relations ( |3^ ) we already 
know that Eq. (|6J) is true Wq + 2i<n — 3 with n = 1. 
For q + 2i = n — 2 Eq. (p3) becomes 



K^'''^~'''\R,0)=j,+,iR)K 

n-4-2i>0 
k=0 



^"''"'*^(i?,0) (65) 

{r,o)k['\r,o) . 



Note that for 2i = n — 2 {n even) one has to resort to Eq. 
( |62| ) and that for 2i = n — 3 (n odd) the sum in Eq. ( |65| ) 
is empty yielding directly 



Kr^ {R, 0) = feiL (R) oc i?C^+2-C"-2 ^ (66) 
/i:W3(ii',0) =7^(i?)fe^(i?) oc i?C"+2-C.-i-3+25„,i ^ 



respectively, which are indeed of the form of Eq. (|64|). 
Now for n — 4 — 2i > the only object in Eq. ( |65| ) whose 
scaling is not known by Eq. ( |64[ ) for g + 2i < n — 3 is 



that corresonding to fc = 0, iiT, 



7r-4-2i) 

i+1 



{R, 0). Hence we 



prove also by induction, but now on i, that Eq.(|64D holds 
for q + 2i = n — 2. Assuming that k\^'^^^^^^ (i?, 0) scales 



according to Eg. (64) we show that it is also true for 
i — 1. By Eq. (|66| ) we already know that this is true for 
2i = n — 4 (n even) and 2i = n — 3 (n odd). We first show 
that for n — 4 — 2i > the term in fc = in Eq. ( |65| ) is 
leading with respect to the terms in higher k appearing 
in this equation. By hypothesis eq. ( |65| ) yields 

(ri-4-2i>0 \ 
_^C2i + 5+C7i + l _j_ y ^ ^C7i + 2-fc+C2i + 4+fc j 



Hence one has to show that for fc > 



Cn+2 ~ Cn+2-fc < C2i+4+fe ~ C2i+4 , 



(68) 



which follows from the Holder inequalities since fc < 
n — 4 — 2i. Secondly, we notice from Eq. ( |67| ) that the 
term in front of the sum scales exactly as a term in fc = 1 
so that by the preceding argument it is subleading with 
respect to the term in fc = for n — 4 — 2i > 0. It follows 
that the term in fc = is the leading one in Eq. ( |67| ) so 
that coming back to Eq. (|63) one has 



K 



(n-2-2i) , 



(R, 0) « k[1-^-^''^ (R, 0)K\''> (i?, 0) (69) 

oc i^?"+2^C2i+2-n+2i+2i5i,_i 



(0), 



Hence Eq. (|^ is true for q + 2i = n — 2 with g > and 
by induction on i for any i > —I. Now by induction on 
n it is also true for any n. We have therefore proven Eq. 
(H) and henceforth Eq. (H). As for Eq. df) it results 
directly from Eqs. (|66|) and (|6^) 



V. THE MULTI-SCALING REPRESENTATION 

The multi-scaling representation of ^2 {R, t) can be 
written as iSl: 



// f}\ 2h+Z(h) / T \ 



where U is the characteristic magnitude of the velocity 
difference across the outer scale of turbulence, is a 
function of the scaled time variable only, and 



(71) 
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The function Z{h) is related to the scahng exponents C„ 
of the nth order structure functions through the saddle 
point requirement 



Cn 



mm\nh 

h 



Z{h)] 



(72) 



To find the scaling exponents associated with the time 
derivatives of ^2iR,T) at r = one computes the n-th 
order partial time derivative of Eq. ( [70| ) to obtain 



in). 



2 (R.r) 

1/2+71 



(73) 



Il-^(,2+n)h+Zih) ^^^^ 



TRJi 



where 



(74) 



At T = this gives 



-r-(") 
2 



(i?,0) 



j.(n) 
12 



(0) 



1/2+11 



Jl^ {2+n)h+Z{h) 



(75) 



Computing the integral at the saddle point in the limit 
R/L ^ and using ([72[) we find 



(76) 



in correspondence with Eq. (p3D. We thus see that the 
continued fraction representation generates the same pre- 
dictions regarding the multiplicity of time scales char- 
acterizing the time correlation functions as the multi- 
scaling representation. We take this as an independent 
evidence for the correctness of the latter. 



VI. CONCLUSION 

In conclusion, we showed that the formally exact con- 
tinued fraction representation of the time correlation 
functions of Belinchier-L'vov velocity differences has the 
same Taylor expansion as the multi-scaling representa- 
tion, at least in terms of the leading scaling exponents 
order by order. 

It should be noted that the continued fraction repre- 
sentation can be used as an approximant for the correla- 
tion function when analytic forms of the time-correlation 



functions are needed. In the lowest approximation one 
takes in Eq. ( |30| ) fci(i?) = 0, producing an exponential 
decay of the correlation function, with a typical decay 
rate 70 (i?). In every successive approximation (^2 = 0, 
fca = 0, etc.) one introduces more and more character- 
istic scales, each one characterized by a different "dy- 
namical exponent" , taking progressively more informa- 
tion about the statistics of higher order correlation func- 
tions into account. 
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